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Thermal fluctuations of the order parameter in an ultrathin triangular shaped superconducting 
structure are studied near Tc, in zero applied field. We find that the order parameter is prone to 
. much larger fiuctuations in the corners of the structure as compared to its interior. This geometry- 



(N 

Oh: 



C/3 



> 

o 



X 



induced localization of thermal fluctuations is attributed to the fact that condensate confinement in 
the corners is characterised by a lower effective dimensionality, which favors stronger fluctuations. 

PACS numbers: 74.25.Bt, 74.25. Uv, 74.78.Na 



I. INTRODUCTION 



a ■ 

O ' Recent progress in growth technology has made it possible to fabricate ultrathin films, consisting of just few 
O , monoatomic layers. Very recently^ superconducting properties of single-atomic layer films made of Pb were reported, 

■ including even the observation of Abrikosov vortices. It is also possible to make nano- and microstructures based on 
such thin films. These nanostructurcs arc generally of much higher quality than those, used in earlier experiments^, 
which, in particular, were not completely free of pinning centers. 

It is well known that the lower the dimensionality of the system the stronger the thermal fluctuations. Ultrathin 

■ superconducting films as well as nanostructurcs based on such films are obvious candidates for the observation of 
^ ' thermal fluctuations. Indeed, it was reported very recently that strong current-induced thermal fluctuations of the 

I , order parameter were observed in superconducting nanowires fabricated in a meander of NbN^. Structures of this kind 
■ O ' are used in highly-sensitive photodctectors. It turns out that the fluctuations arc a main source of dark-count events 
C , in such photodctectors: this means that their effect is actually parasitic. Fluctuation phenomena were interpreted in 
^ ' terms of a thermally-activated entry of vortices, as well as by unbinding of vortex-antivortex pairs^. In another recent 
experiment^, no hysteresis for vortex penetration and expulsion was observed in the case of a nanoisland of Pb, which 
was so small that it could accommodate only one vortex. In a similar experiment^, performed at lower temperature, 
some hysteresis was detected, but the width of the hysteresis region was significantly smaller than expected from 
theory. A theoretical explanation for these phenomena was very recently suggested by one of us in terms of a thermal 
suppression of the surface barrier for vortex entry/exit, which might occur in superconducting nanoislands made of 
Pb in the regime of the ultimate single vortex confinement^. Previously, thermal activation of vortices over the surface 
' barrier was demonstrated to be possible in high- Tc low-dimensional structures^. 
[ Geometry itself can play an important role in fluctuation phenomena. The shape of islands, prepared by the method 
. of evaporation in ultrahigh vacuum, depends strongly on their sizes. They grow according to the Stranski-Krastanov 
' scenario, i.e., starting from well- facetted nuclei^. As a result, islands with lateral dimensions of nearly 100 nm and 
y—i . smaller have a hexagonal shape, while larger islands tend to have a triangular shape [see, e.g.. Fig. 1(a) in Ref4i]. At 
ILJ ■ the same time, it is known that the triangular geometry can lead to very peculiar consequences for the superconducting 
condensate confined in this geometry. For instance, it was shown in Refsi^i^ that stable vortex-antivortex molecules can 
nucleate, when a homogeneous magnetic field is applied. The aim of the present paper is to explore how triangularity. 



JH ] or more generally the presence of sharp corners in a thin superconducting structure, influences thermal fluctuations 
- - < of the order parameter. The motivation to study the triangular geometry, apart from possible links with experiments, 
is that one can expect thermal fluctuations to be stronger in the corners of a triangular structure as compared to its 
interior, since the superconducting condensate in the corners is strongly conflned by samples's borders. This implies 
that the system somehow is characterized by lower effective dimensionality than in the interior. It is however not 
evident a priori if this rather general argument leads to any noticeable effect for the fluctuations of the order parameter 
in a triangle, so that more careful inspection is certainly desirable. In order to reveal the effect of geometry on spatial 
localization of thermal fluctuations, we restrict ourselves to the case of zero applied fleld, but we consider temperatures 
both below and above Tc- We apply the method of small- amplitude oscillations within the Ginzburg-Landau theory, 
and we do not consider the possible excitation of vortices. 

The paper is organized as follows. In Section II we formulate our model. In Section III we present our results for 
several correlation functions both below and above T^. We conclude in Section IV. 
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II. MODEL 

The derivation presented in this Section appUes for the case of temperatures below T(.. It is straightforward to 
adopt it for the case of temperatures above T^. 

A. General formulation 

We start with the dimensionless Ginzburg-Landau functional for the superconducting energy of the island of thick- 
ness d, in zero applied field: 



F = ^^£,{Tfd I d\ (- l/l^ + \ \.ft + |V/r ) , (1) 



where integration is performed over the cross section of the nanostructure, / is the dimensionless order parameter. 
All distances are measured in units of the temperature-dependent coherence length £,{T), and Bc(T) is the thermo- 
dynamical critical field given by 

27rV2e(r)A(T)' 

We consider ultrathin islands, with d C(T'), so that the problem is two-dimensional. The boundary condition for 
the order parameter at each of the three edges of a sample is taken in its usual form for the case of a superconduc- 
tor/vacuum interface 

If = 0, (3) 
an 

where n is the unit vector in the direction perpendicular to the edge. 

Let us now estimate the ratio G{T) of the energy needed to suppress the order parameter to zero within the volume 
~ ^{T)^d and the thermal energy /cbT, using realistic values for all parameters for typical ultrathin nanostructures 
made of Pb that have been realized in recent experiments. We define G{T) as 

2fiakBT 



For the penetration depth, we use the usual expressioi*i^ for dirty superconductors, A(T) ~ 0.62Ao y where 
Ao — 40 nm is the penetration depth in bulk Pb, Z ~ 2(i is the quasiparticle mean free part, and Tc = 7.2 K. Then, 
for a film of thickness 2 nm (similar to Ref.— and nearly 3 times larger than in Ref.— ), operated at T = 0.75Tc, we 
estimate G(T) - 10. 

B. Fluctuative modes 

Next, we expand the order parameter / around its mean value in the absence of fluctuations 

f = l + 6f. (5) 

After substituting Eq. (5) into Eq. (1), we obtain the expression for the increase of superconducting energy due to 
thermal fluctuations 



SF = -^j-^aTYd J d'r i^\Sfr + -iSr + Sr') + |V5/r j , (6) 

where we have kept only terms up to quadratic ones in 6f and Sf*. 

Obviously, Sf can be represented as a superposition of plane waves. However, Sf for each fluctuative mode must 
satisfy the boundary condition given by Eq. (3), which is a very strong restriction imposed by the geometry of our 
problem. In order to circumvent this difficulty, wc will use, with certain modifications, an approach applied very 
recently in RefJ^ for the triangular dot of graphcne, for which exact electronic wave functions were found. The major 
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FIG. 1: Schematic image of a triangular structure. Numbers depict points, the correlation functions being calculated along 
the straight lines connecting these points [see in the text] . 



difference with the graphene dot problem is in the boundary condition, which is vanishing of the wave function at the 
triangle's border, instead of the vanishing of its first derivative [Eq. (3)]. Following Ref.— , as a first step we focus on 
the sector, confined between edges I and II [see Fig. 1] and we find the form of 5f, which satisfies boundary condition 



[Eq. (3)] at these two edges. Consider a plane wave with the wave vector ki 

ipi = exp(— ikir). 



(7) 



ky). 



It is then easy to see that the sum 



When reflected from edge I, it is converted into wave ip2 with k2 = {kx 

of -01 and 02 satisfies Eq. (3) at the edge I. This is in contrast with Reffii, where the difference between tpi and tl'2 
was taken due to the different boundary condition. When these two waves are reflected from edge II, two new waves 
appear with ks = —^{kx + \/3ky,—\/3kx + ky) and kg = —^{kx — V3ky,—\/3kx — ky). After the reflection from edge I, 

these two ones give rise to two more waves with ks = —^{kx — VSky, y/Skx + ky) and k4 = —■^{kx + V^ky, VSkx — ky). 
Finally, after reflecting from edge I, these last waves do not lead to any new waves. Then, the sum of the six wave 
functions ipi + ...+ ipe satisfies the boundary condition, given by Eq. (3), both at edges I and II, whatever ki is. The 



same boundary condition at edge III leads to quantization rules for {kx 
It is easy to see that the allowed {kx,ky) split into two branches: 



ky), which differ from those found in Refi^ 
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where L is the dimensionless length of the triangle's edge, while m and n are arbitrary integer numbers. 
For these two branches, eigen energies are given by 
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Wave functions corresponding to different (m, n) but the same brancfi are not necessarily unique. By a proper choice 
of m and n, we obtain the following basis of normalized and unique wave functions 

= -L(^W + ... + 4i)), n>\m\, (14) 
Vmi = 1, ?^ = m = 0, (15) 
where ■>p[^\..., ip^^^ depend on {ki^\ky^^), and 

= i= [4'^ + ... + 4')) , n > \m + 1/2| , (16) 

where V'e^^ depend on {k'i'\k!if' ) . 

We can expand the fluctuating contribution to the order parameter in the constructed basis: 

n,m,Q 

where a = 1,2 stand for the two branches. We now substitute expansion (17) in Eq. (6). After performing the 
integration over the triangle's area, we arrive at the following expression for the statistical sum 



Z^W I d{Re{ctl))d{Im{c^;^l)) exp | -G{T)^ [(i?e(c(;^)))2(i?(^) + 2) + (/m(cM ))'i?(?, 



It contains a product of Gaussian integrals, which can be easily evaluated analytically as 



Z=\[ (18) 

lJeI^I{e'^^1 + 2)G(r)V3 



C. Correlation functions 



In order to study the spatial localization of fluctuations in the triangle below Tc, we analyze separately phase and 
density fluctuations, since they behave in different ways. Namely, these are fluctuations of the phase of the order 
parameter which are responsible for the loss of the long-range order—. 

It is straightforward to see, from Eq. (5), that the phase x of order parameter can be expressed as 



The quantity we are interested in is the correlation function 

A>.(ri,r2)==([x(ri)-x(r2)]')^ 



(19) 



(20) 



In an infinite isotropic system this quantity depends only on the distance |r2 — ri|. It diverges logarithmically at large 
distances for two-dimensional systems, expressing the loss of long-range order. However for a confined system with 
anisotropic geometry, Kph{Yi,T2) depends on both ri and r2 — ri. 

Using Eqs. (17) and (19), we arrive at the following expression for Kph{Ti^Y2) 



K,n[v,,V2) = \Y. ((/m(4"2))^)^|i?e 

m,n,a ^ 



(ri,r2) 



(21) 



where 



(22) 
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Thermally-averaged {Re{cmn)Y {Im{cmn)Y can be found analytically as 

2 

G(r)L2V3 + 2 



(i?e(4"^))')^ = ^, (23) 



{Im{c[ZL)f) = V r r (24) 

By substituting Eqs. (23) and (24) into Eq. (21) and performing the summations, we calculate Kph{Yi,Y2) along 
various lines inside the triangle. This sum, however, is divergent at large m and n, and this divergency has to be cut 
in a standard way [see, e.g., p. 336 of Ref.— ] at wave vectors k ~ 1/^(T). We found that the dependence of the final 
results on the particular choice of the cutoff value is rather weak, within 10 %. 

In a similar way, we can express the fluctuation part 811^ of the density of superconducting electrons |/|^ in terms 
of 5f and 5f* 

5n, = 5f + 5f*. (25) 

By using the developed approach, it is also possible to find the density-density correlation function, defined as 

A'den(ri,r2) = {5ns{Yi)5ns{v2))rp. (26) 

Thermal fiuctuations of the order parameter also occur at temperatures above T^. In this case, the equilibrium 
order parameter is equal to zero, which leads to changes in Eqs. (5), (18) and therefore also in Eqs. (23) and (24). 
We omit the derivation and we will present our final results in Section III for the correlation function for the order 
parameter itself, without separating the density and the phase 

A-(ri,r2) = (/(ri)*/(r2))T- (27) 



Above Tc, it is also convenient to normalize all the distances by ^(T), where ^(T) = ^Jh? /2m |a(r)|. 



III. RESULTS AND DISCUSSION 



In this Section, we discuss the behavior of the correlation functions in the triangle both below and above Tc. 



A. Temperatures below Tc 



Our results for the dependence of the phase correlation function Kph{Yi,Y2) on \y2 — Yi\ along different paths are 
shown in Fig. 2 for a triangle with L = 8 [several hundreds of nanometers for the case of a nanostructure made of Pb 
at T = 0.75Tc, G{T = 0.75Tc) = 10]. Different curves show the behavior of Kph{Yi,Y2) along straight lines connecting 
various points inside the triangle, which are depicted in Fig. 1. Curve 1 — > 2 corresponds to the path between the 
two corners of the structure. Curve 1^3 gives Kph{Yi,Y2) along the bisector starting from the corner. Curves 4 — )• 3 
and 5 — >■ 3 correspond to lines in the interior region of the structure [lengths of segments 1 — 4 and 2 — 3 are L/2, 
while lengths of segments 4 — 5 and 5 — 3 are the same]. Curve 6 — >■ 7 shows Kpii{Yi,Y2) across the bisector [lengths 
of segments 1 — 6 and 1 — 7 are L/5]. 

The shape of all these curves is similar and it generally resembles the behavior of the same correlation function for 
an isotropic infinite two-dimensional system. Namely, it first grows rapidly with increasing \y2 — ri| and then starts 
to demonstrate a smoother behavior. There are, however, obvious differences between the curves, both quantitative 
and qualitative. First, it is seen that the initial increase of A'p/i(ri,r2) is much larger [several times] for those paths 
which originate from the corner [curves 1—^2 and 1—5-3], where, according to our initial guess, fluctuations are 
stronger. Thus, we can conclude that this guess is justified. Second, Kph{Yi, Y2) along the curve 1 — 2 starts to rapidly 
increase again, when approaching the other corner, while the such an increase for other curves is not so pronounced. 
This feature implies that the order parameter between different corners is more decoherent, which again supports 
our hypothesis. For the range of parameters used here, we see that the maximum Kph{Yi,Y2) is nearly 0.4, which 
implies that the coherence between the order parameter inside different corners starts to be lost. 

Figure 3 gives the density-density correlation function KdeniYi, Y2) for three directions, along paths 1 — s- 2, 4 — )• 3, 
and 5 — > 3. We again see that the correlation function in the corner [1 — )■ 2] is much larger than the same function 
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FIG. 2: (Color online) Dependence of the correlation function _K'ph(ri,r2) on [r2 — ri[ along different paths in the triangular 
structure below Tc- Red lines show Kph{v\,V2) along the paths, which originate from the corner. Blue lines correspond to 
paths in the interior of the structure. Green dashed line corresponds to the path across the bisector. 
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FIG. 3: (Color online) Dependence of the correlation function Kden{T^i,V2) on |r2 — ri| along different paths within the 
triangular structure below Tc- Red line shows Kdmijci, r2) along the path, which originates from the corner. Blue line 
corresponds to path, which starts from the edge towards the interior of the structure. Green line shows Kdm{ri,Y2) starting 
from the interior to the edge. 



in the interior [4—5-3 and 5 — > 3]. This shows that not only the phase of the order parameter, but also the density 
of superconducting electrons fluctuates stronger in the corners. It is also of interest to note that the comparison of 
Kden{Yi,Y2) aloug paths 4 — >■ 3 and 5 — 3 shows that fluctuations are stronger near the triangle's surface as compared 
to the central part of the nanostructure. 

The strength of thermal fluctuations at the given point ri may be characterized by the density-density correlation 
function i^de„(ri,r2) with |r2 — ri| ^ ^(T) averaged over all such values of rjJ^ii^. Although i^Tden (r i, is usually 
peaked within nearly 1.5^(T) in the corners for various triangle's sizes, such an averaging implies that the typical 
width of the geometry-induced fluctuating region in each of the three corners is broader: we estimate it to be several 

In order to visualize the effect of geometry on thermal fluctuations, we construct a "phase diagram" for the case 



7 



0.6 
0.5 

T- 0.3 

0.2 
0.1 

1 2 3 4 5 

distance [f (0)] 



FIG. 4: "Phase diagram" of the triangular structure below Tc indicating the region where fluctuations are more pronounced. 

of a triangular Pb nanostructure with the edge length 15^(0) and thickness d = 5 nm. We take 1^-3 path [along 
the biseetor] and for eaeh point r at this path wc calculate •\//^deri(r, r), which shows the average fluctuation of the 
density of superconducting electrons. We have chosen the following qualitative criterion: if this quantity exceeds 0.15 
at a certain temperature, we enter the fluctuative regime. Our results for the crossover between low-temperature and 
fluctuative regimes are presented in Fig. 4 in the plane of the distance along the 1 — >■ 3 path and 1 — T/Tc- We 
indeed see that the characteristic temperature for the crossover in the corners strongly depends on the position. Of 
course, the transition between the two domains in Fig. 4 is not abrupt. It is more appropriate to talk of different 
regimes, which are separated by a crossover region. Thus, the suggested picture is very qualitative and it is given for 
illustrative purposes to supplement an idea we want to convey. 

B. Temperatures above Tc 

The order parameter correlation function A'(ri, above Tc, defined in Eq. (27), has to decrease^^ with increasing 
|r2 — ri|. The results of our calculations for A'(ri,r2) along different paths are presented in Fig. 5, where the 
correlation function is plotted for T = 1.2hTc for the triangle with length L = 8 [in terms of ^(T)]. K{ti,Y2) is 
normalized by its value A'(r,r) in the point 5, which is located in the central part of the nanostructure ["bulk"]. 
We indeed see the expected rapid decay of Ar(ri,r2) with increasing |r2 — ri|, but again the correlation function 
calculated along the path, which originates from the corner [1—^2], is rather different from those, which correspond 
to lines inside the interior of the structure. Namely, the average value of the density of superconducting electrons 
induced by thermal fluctuations is much larger in the corners. Fluctuations near the surface are stronger than far 
from the surface [curves 4 — > 3 and 5 — 3]. The width of the geometry-enhanced fluctuating region in the corners is 
again estimated to be several i{T). 

The comparison of our results for the correlation functions below and above Tc shows that superconductivity is 
getting suppressed, below Tc, starting from the corners, while, above Tc, the order parameter preferentially nucleates 
again in the corners. However, there is no contradiction, since if we consider a triangle with fixed sizes and if we start 
to increase T towards Tc, ^(T) becomes infinitely large at Tc. This means that, close enough to Tc on both sides of 
the transition, the triangle is actually in the zero-dimensional regime, so that the order parameter does not vary in 
space. In other words, the fluctuating region in the vicinity of Tc expands from the corners to the whole structure, 
while by tuning temperature away from Tc, one can make thermal fluctuations stronger in the corners than in the 
interior. Of course, the length of the triangle's edge should be much larger than ^(0) in order to see such a spatial 
localization of thermal fluctuations. 

The geometry- induced localization of thermal fluctuations is expected to occur in nanostructures of various shapes, 
not only triangular ones. The general tendency is obvious: the sharper the corner, the stronger the fluctuations in 
this corner. For instance, in the corners of square-shaped samples, the width of the fluctuative region can be expected 
to be much smaller than in the corners of the triangular structures. 




fluctuative regime 
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FIG. 5: (Color online) Dependence of the correlation function if (ri, r2) on \v2 — ri | along different paths within the triangular 
structure above Tc. Red line shows K{r\,Y2) along the path, which originates from the corner. Blue line corresponds to path, 
which starts from the edge towards the interior of the structure. Green line shows _K'(ri,r2) starting from the interior to the 
edge. 



IV. CONCLUSION 

We have studied thermal fluctuations of the order parameter in quasi-two-dimensional superconducting structures 
of triangular shape. We considered the case of zero applied field both below and above Tc- It was shown that the order 
parameter exhibits much larger fluctuations in the corners of the structure and the width of such a geometry-induced 
fluctuating region is several ^(T). This unusual behavior can be attributed to the fact that the confinement of the 
condensate in the corners lowers locally the effective dimensionality of the system, thus making thermal fluctuations 
more pronounced. The condensate in the corners can serve as a source of thermal noise in superconducting devices 
based on ultrathin nano- and microstructures. 
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